Artifacts of opinion dynamics at one dimension 
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The dynamics of a one dimensional Ising spin system is investigated using three families of local 
update rules, the Galam majority rules, Glauber inflow influences and Sznadj outflow drives. Given 
an initial density p of up spins the probability to reach a final state with all spins up is calculated 
exactly for each choice. The various formulas are compared to a series of previous calculations 
obtained analytically using the Kirkwood approximation. They turn out to be identical. The 
apparent discrepancy with the Galam unifying frame is addressed. The difference in the results seems 
to stem directly from the implementation of the local update rule used to perform the associated 
numerical simulations. The findings lead to view the non stepwise exit probability as an artifact 
of the one dimensional finite size system with fixed spins. The suitability and the significance to 
perform numerical simulations to model social behavior without solid constraints is discussed and 
the question of what it means to have a mean field result in this context is addressed. 



In this Letter, we reinvestigate the dynamics of an ex- 
tremely simple opinion dynamics model, a slightly mod- 
ified version of the original Sznadj outflow rule [l| pro- 
posed by Slanina et al which was shown to exhibit 
a rich variety of behaviors. In particular, when calculat- 
ing analytically the probability to reach a uniform final 
state, they found it is a continuous function of the initial 
magnetization. The result was confirmed by Monte Carlo 
simulations. Starting from a generalized voter model an 
independent and simultaneous paper by Lambiotte and 
Redner Q recovered both results, analytically and nu- 
merically. The existence of a continuous exit probability 
contradicts the prediction of the Galam unifying frame 
(GUF), also denoted as the general sequential probabilis- 
tic frame (GSPF) Q, which yields a threshold stepwise 
function for the exit probability. 

We show that these analytical results are obtained at 
once by treating one single triplet instead of the total sys- 
tem. Alternatively, a straightforward one shot extended 
voter model is found to reproduce the same results. The 
two treatments are a direct application of classical mean 
field technics, which indeed are based on a "one-site ap- 
proach" [Bj. In addition, a simple sequential update is 
found to be sufficient to recover the results via a Monte 
Carlo simulation. 

Based on the results of this mean-field approach ap- 
plied to triplets, we arrive at three conclusions: (i) The 
Kirkwood approximation Q treatment of the Slanina 
chain (and, by extension, and the Lambiotte and Redner 
treatment) reduces to a simple mean field approach, (ii) 
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Contrary to several author claims, the GUF does go be- 
yond the mean field. Indeed a one step GUF restricted 
to an average single interactive group yields the mean 
field level, (iii) The Monte Carlo simulations, which were 
found in perfect agreement with the analytical formula, 
must have been biased either by finite size effects or by 
using an inappropriate sequential update. 

To set the stage for the demonstration, we recall that, 
in statistical physics, microscopic interactions are given 
by defining the Hamiltonian of the system and the as- 
sociated macroscopic properties obtained from analyti- 
cal calculations or from simulations that are expected to 
be independent of the choice of the dynamics used to 
implement those interactions. Moreover, while perform- 
ing simulations, a special attention is given in the actual 
choice of the update procedure, either sequential or si- 
multaneous or random, to avoid results which could be 
artifacts of the procedure used, with the system being 
trapped in some specific configurations. 

However, in opinion dynamics, most models focus on 
the choice of the ripdate rule without formulating the ac- 
tual interactions [3] . Some exceptions exist [1, Q • It is 
similar in economy with game theory where only the out- 
comes are give n without an explicit formulation of the in- 
teractions [lOl • The social content of the model is exhib- 
ited in the setting of the update rule, which is supposed to 
incorporate the complex and unknown mechanisms of the 
underlying cognitive interactions. Understanding what a 
mean-field approach is under these circumstances can be 
difficult. This is probably the cause for the confusion 
about above approaches, which turn out to be approxi- 
mations. 

The fact is that a great deal of models have been pro- 
posed to describe opinion forming, each one displaying a 
specific update rule [11, [TT| - [l6| . Usually, this makes it dif- 
ficult to perform analytic calculations and most of the re- 
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suits are obtained from numerical simulations where the 
rule can be implemented. The significance of the model is 
then discussed on the sole basis of the numerical results, 
often obtained from small size samples, usually of the 
order of 100 or 1000, rarely more. Such a practice, com- 
bined with the almost impossible making of real experi- 
ments, renders it difficult to assess the social relevance of 
a model without solid ground. Even the validity of the 
numerical results cannot be certain since they are not 
constrained by some general accepted principles, which 
have to be obeyed like in physics. Another difference 
is that in statistical physics interactions are simultane- 
ously activated among the spins although they usually 
involve only nearest neighbors while in social systems in- 
teractions are restricted to separate groups of agents at 
a time via local interactions among the agents. 

Here, we first outline the basic features of the modi- 
fied Szandj model with the approximations used to solve 
it. The corresponding results obtained by analytical cal- 
culations and numerical simulations are listed. Then we 
present Galam exact derivation of the same results and 
seek to understand three basic issues: (i) What is the 
physical implications of that identity? (ii) What is the 
social significance of the associated Monte Carlo simu- 
lations, (iii) What is the validity of an exit probability 
which is not stepwise? We demonstrate that, while GUF 
provides a mean field result if iterated only once, it is 
not mean field when iterated more than once and to sev- 
eral local groups. Applied to the Slanina model the one 
step averaged GUF reproduces the same result reported 
previously 0, [1]. 

The Slanina et al model and results. Consider N agents 
on a linear chain at sites denoted respectively 1, 2, N . 
A two state variable Si = ±1 with i = 1,2,...,N is at- 
tached to the agent at site i, which represents its current 
opinion, cither -1-1 or —1. Given a distribution of individ- 
ual opinions the system evolves by repeated applications 
of the following procedure: (a) Two neighboring agents 
are selected randomly, (b) In case they have different 
opinions, nothing happens and a new pair is selected, (c) 
If they do share the same opinion, one neighbor of the 
pair is chosen randomly with equiprobability. (c) The se- 
lected neighbor adopts the common opinion, if it didn't 
already agree to that, (d) The process is repeated till all 
agents hold the same opinion. While two opinions can 
flip simultaneously in the original Sznadj formulation, 
here one spin is flipped at most in one step. 

Given a distribution of -1-1 and —1, the probability P_|_ 
to have all agents sharing the opinion -t-1 is calculated 
as a function of the density pq of -1-1, which are initially 
present in the chain. It is called the exit probability. 
The number of updates required to reach the state with 
all agents aligned is also calculated. 

The calculation of the exit probability is carried on ap- 
plying a decoupling approximation to truncate the equa- 
tion hierarchy for multi-point spin correlations. It is used 
in various contexts with different names. The one con- 
sidered was introduced by Kirkwood to study the statis- 
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FIG. 1: The variations of botli P+ from Eq. U]) and pi from 
Eq. (??) are exhibited as a function of po- 



tical mechanics of fluid mixtures [6| . It is similar to the 
Hartree-Fock approximation but is more limited since it 
cannot be systematically improved with diagrammatic 
techniques. A second approximation is applied with the 
assumption of low decay of correlations. Nearest neigh- 
bor correlations are extended up to next next nearest 
neighbors. After elaborated manipulations the following 
exit probability is obtained with 



- 2po + 1 



(1) 



The function is shown in Figure © and was found in 
perfect agreement with Monte Carlo simulations of chains 
of sizes up to 10'^ agents with lO** sampling Identical 
results were obtained in Q with 100 agents averaged over 
5 X 10"^ realizations. 

The novel treatment. In order to get a mean field result 
to the model, we first reformulate it by treating a triplet 
instead of a pair plus a random neighbor. Pick at random 
one site Si. Then complete the triplet with equal prob- 
ability, either on its left with Si^2iSi-i or on its right 
with Si+i,Si+2- In the first case, if 5^-2 = Si-i = a, 
then Si = a and in the second case if 5'^+ 1 = S',;+2 = b, 
then Si = b with a,b ~ + or -. Otherwise for 5^-2 ^ Si-i 
or Si+i Si+2, nothing happens and a new site is se- 
lected. We consider that one interaction has happened 
only when one site does change its value. 

Accordingly, the probability to have a selected site i to 
end is just + ^p^ = Pq, the probability to have a 
pair of + on the right or on the left. By symmetry the 
probability to have a selected site i to end -, is just (1 — 
Po)^, the probability to have a pair of -. However, since 
the case of a mixed pair Si^2 = —Si-i or 5*^+1 = —Si+2 
is dismissed with a new site Si being selected, the sum 
of above two probabilities is not equal to one. Therefore 
they must be rcscalcd by their sum to get normalized to 
one. Given above update rule the normalized probability 
to have a site updated to -|- is thus. 
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The same result is obtained by picking randomly a pair 
of adjacent sites. If they are in different states, select 
another pair. The probability of this event is 2p{l ~p). 
Otherwise, if they share the same state, i.e. + or - with 
respective probabilities p^ and (1 — p)^, select the left or 
right adjacent site and align it along the pair common 
state. The normalized probability to get an updated + 
site is thus given also by Eq. ([5]). 

At this stage, it is worth stressing that Eq. ([2]) is 
identical to Eq. (jlj . However while P+ is the probability 
to end up with all sites at + after a series of repeated 
updates, pi is the probability to get one site + after one 
update for one selected site and triplet. The next step is 
then to calculate the exit probability which results from 
successive iterations of Eq. ([2]) for all sites. 

One possible way to implement the calculation is by 
using the following sequential update procedure: Pick up 
one pair of adjacent sites, if they arc heterogeneous, pick 
another pair. Once the selected pair is homogenous, align 
a random adjacent site along the common state, either 
on the left or on the right. To illustrate the procedure, 
consider the case where the pair is ++ with 

...•••(++)••••...; (3) 

Select the adjacent site to be updated on the right (upper 
line) or on the left (lower line), 

( ...•••(+ + •)•••... 

-•••(++)••••■■• ^ < ; (4) 

[ ...••(• + +)••••... 

Update for instance the right upper line). It yields, 

...•••(+ + +)•••...; (5) 

Move the triplet one site left or right as follows, 

( ...•••+(+ + •)••... 
...•••(+ + +)•••... ^ <^ ; (6) 

[ ...••(• + +) + •••... 

Update one of the triplet with one of the two chains, 

( +(+ + +) • «... 

{ ; (7) 

[ ... ••(+ + +) + •• «... 

Repeat the process by shifting one site left or right the 
triplet as above. The whole chain will end up with all 
sites in state +. Accordingly, the only random step is 
the selection of the initial pair, then the process is de- 
terministic with respect to the final alignment. It proves 
that the sequential update procedure leads to P-f ~ pi . 

The elaborated Kirkwood decoupling approximation 
to truncate the equation hierarchy for multi-point spin 
correlations is thus identical to above straightforwards 
calculation. 

The one step mean field voter model extension. A sim- 
ilar analysis can be performed by using the voter model 
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FIG. 2: The successive iterations of Eq. (??) with 
pi,P2,P3, ...ps as a function of pt- The drive towards a step 
function for the exit probability is clearly seen. 



as basis. Voter model is based on a conservative dynam- 
ics using a pair directed interaction (T^]. Select one site 
Si at random, pick a second one Sj again at random, 
then update the first site putting Si = Sj with Sj un- 
changed. Within the GUF, the voter model was shown 
to be the frontier between the ordered phase, charac- 
terized by the coexistence of a majority and a minority, 
including the extreme limit of a zero minority, and the 
disordered phase where a perfect equilibrium between the 
two opinions prevails 

A mean field treatment of the voter model @ consists 
of treating one site exactly while the interacting sites are 
taken at their mean value. As the system is at zero tem- 
perature, there is no thermal variations in the sites that 
are treated exactly and we need to treat the problem dif- 
ferently. With no probabilistic distribution of spins due 
to thermal noise, the only probabilistic aspect left is that 
of random distribution of spins on the sites. The only 
two possible equilibria are when the pair site-external 
field (external site) have the same spin. As both have 
the same probability of being up, initially, the chain is 
all + with the probability Pq and all - with probability 
(1 — po)^. However since the case of a mixed pair is dis- 
missed, both probabilities should be normalized, which 
reproduces P_|_ = pi given by Eqs. ([1] This restric- 
tion to one step update of one site produces a mean field 
result. Our demonstration enlightens the mean field na- 
ture of the previous derivation. 

The only difference between the various calculations 
is the number of Monte Carlo steps needed to reach the 
final sate with all sites sharing the same opinion. These 
numbers are a direct function of the procedure chosen to 
implement the updates. They are not robust but could 
be used to make social predictions with respect to setting 
some procedural frame for decision making. Indeed, it is 
an efficient way to gain a lot of time to reach the same 
final decision. 

Applying the GUF. Above results help to understand 
the discrepancy noticed by Slanina about the GUF, 
which is expected to yield a step function for P^, while 
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other approaches have provided a continuous function. 
We have seen that the results of those approaches can 
be obtained from the analysis of a very small number of 
interacting sites, by using the GUF approach and can be 
seen, in each case, as a consequence of the method used to 
obtain the result, instead of the correct thermodynamic 
solution. 

In the case of the modified Sznajd rules, the Kirkwood 
approximation produced the same result for P-^ as the 
probability pi that the first spin to change in the system 
will change to +. And we have shown that in a simu- 
lation, depending on the update procedure, that proba- 
bility can be the actual final probability that the whole 
system will end as +, but that this was obtained as an 
artifact of the implementation rule and not a correct re- 
sult. Another possibility is that the observed simulation 
results could be a consequence of finite size problems. 
Indeed, the exact result for TV = 3 sites is very easy to 
obtain as the system is completelly deterministic. It is 
basically the probability that the randomly drawn spins 
will have a majority of -|-s and we get P+ = — 2pq. 
This curve is further away from the step function than 
Kirkwood approximation. As the N becomes larger, P_|_ 
should get closer to the Kirkwood solution in its path to 
a step function and this might be another cause for the 
observed behavior of the simulations. 

In summary, we have demonstrated that Slanina and 
Lambiotte results can be obtained from a mean-field ap- 



proach. On the other hand, the GUF treatment of the 
same problem does provide a continuous function, qual- 
itatively similar to the ones obtained by the Kirkwood 
approximation, if one uses only one interaction. How- 
ever, the sites will keep updating their spins and, after a 
few iterations, it is easy to see that the GUF approach 
tends to a step exit probability. The additional subse- 
quent steps drive the system further apart from mean 
field-like result and this shows that GUF is actually not 
a mean field result. This helps understand better the 
subtleties of extending the concept of a mean field result 
to a zero temperature problem. 

We conclude that when the hamiltonian of the system 
is not properly defined and, instead, update rules are pro- 
vided, special care is needed in order to avoid unwanted 
artifacts of the implementation details. Finite size ef- 
fects with the system still out of equilibrium might be 
far more important than simulations might indicate. It 
is true that social systems are not so large that finite size 
systems will be completelly avoided and, in that sense, 
the previous results of Slanina and Lambiotte can have 
significative consequences in real systems. But, from a 
Physics perspective, they are probably consequences of 
updates rules and the finite size effects. 
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